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ABSTRACT 

A  simple  primitive-equation  model  for  numerical 
prediction  of  wind  fields  and  temperature  fields  at  S00  and 
400  mb  is  developed  and  programmed.   Numerical  experiments 
are  conducted  with  this  model  in  an  effort  to  suppress  the 
spurious  instability  encountered  in  solution  of  the 
primitive  equations .   The  model  proves  to  be  very  sensitive 
to  spatial  finite-differencing  schemes.   The  schemes 
developed  in  this  study  tend  to  suppress  the  instability 
resulting  from  truncation  and  round-off  error.   Instability 
resulting  from  unrealistic  boundary  conditions  for  a 
limited  forecast  region  is  reduced  by  not  allowing 
advection  across  the  boundaries  and  appears  to  be  con- 
trolled by  the  finite-differencing  schemes  used  in  this 
study. 

The  writers  wish  to  express  their  appreciation  to 
Professor  George  J.  Haltiner  of  the  U.  S.  Naval 
Postgraduate  School  for  his  guidance,  contributions  and 
encouragement  in  this  work.   Appreciation  is  also  expressed 
to  the  personnel  of  the  U.  S.  Navy  Fleet  Numerical  Weather 
Facility  for  their  cooperation  during  the  preparation  of 
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1.   Introduction. 

It  is  generally  recognized  that  solution  of  the 
primitive -equations  offers  the  most  general  and  logical 
approach  to  numerical  prediction.   Previous  attempts  at 
solution  of  the  primitive  equations  by  numerical  methods 
have  met  with  limited  success.   However,  interest  has 
recently  been  stimulated  in  the  primitive  equations  by 
unpublished  reports  indicating  that  Arakawa  [7J  has 
succeeded  in  developing  a  stable  solution.   The  purpose 
of  this  study  is  to  develop  and  program  a  simple  3- level 
primitive  equation  model  and  perform  some  numerical 
experiments.   The  main  effort  in  these  experiments  is 
devoted  to  the  suppression  of  spurious  instability 
encountered  in  numerical  solution  of  the  primitive 
equations.   It  is  hoped  that  this  instability  can  be 
suppressed  by  a  finite-differencing  scheme  that  will  in 
the  mean  conserve  the  sum  of  potential  and  kinetic  energies 
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2.   The  Model 

Phillips  dlj  describes  a  coordinate  system  where  the 
vertical  coordinate  p  in  the  x,  y,  p,  t-system  is  replaced 
by  an  independent  variable  V,  where  ^7~=  p-Pt.  The  following 
relation  then  holds,  where  q?  can  be  x,  y,  or  t : 

The  subscript  p  denotes  derivatives  taken  along  a  constant 
pressure  surface .   The  variable  V  ranges  monotonically  from 
zero  at  the  top  of  the  atmosphere  to  unity  at  the  ground. 

If  the  above  coordinate  system  is  utilized  the 
horizontal  equation  of  motion  becomes 

£V  +(v-7)v+-  trdy  +•  v$  -.si&fc  vtt  + 

FF  +  £  k'/V  =    O.      (2) 

In  this  study  it  is  assumed  that  the  top  of  the 
atmosphere  P^  is  at  200  mb.  The  atmosphere  is  divided 
vertically  by  five  levels  as  follows: 

^  =0 p=200  mb 

v^  =1/4 p=400  mb 

$1   =l/2_ p=600  mb 

^5  =3/4 p=#00  mb 

%  =1 SFC 
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For  this  model  7/"*,  and  ^T  will  be  approximated  as  follows 

Vertical  velocities  are  identically  zero  at  the  top  and 

•     •  • 

bottom  of  the  atmosphere  (^£=0,^1=0),  reducing  U~  and 


^   to 
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If  the  preceding  simplifications  are  utilized  and 
equation  (2)  is  applied  at  the  ^  and  ^C  levels,  the 
following  equations  ensue 


J2L  M^  7/T  +  /f  ffft**;  -  o  (4) 

77-  aT 


7T  dT 


It  is  assumed  in  (4)  and  (5)  that 

a  ^r  d  ^         a  t  a<t- 
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Moreover  the  term  V" d $|  in  (4)  may  be  replaced  by 


dT 


^"d^r  (*&5I)  -  cb.  iSI 


Since   ^T  =0  and  c*js[7      =1, 


i    ) 


dT 


T" 


and  equation  (4)  becomes 


i^i  +(y.7)x+4?r(v3-^)  +  ?<fct(4>, 


i£V 


Similiarly  at  the  ^  level 

^^       AT" 

and  equation  (5)  becomes 


(6) 


2d<T 


>3 


4 


^d?4-^^]^V7T  +  l|  +  l?X£^  =  o 


.  (7) 


Equations  (6)  and  (7)  may  also  be  expressed  in  component 
form  as  follows : 
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^0 


^-(♦.-Sfcteff+fc 


(d) 


/iV 


(9) 


I  u3  ^  +  ^  4£  -^  v5 1  -  -SL  Kr *0-  4&  - 


iu3   .  _ 


(10) 


av,  - 


pU3  ^13  +  ^^3+  T^-fcj  "I  -JGtC^  - 


*)- 


cfc&3 

<>3 


(11) 


The  equation  of  continuity  expressed  in  V"  coordinates 
becomes 


<UT+  V-(7TV)  +  j^iTTT)  ^O  • 


(12) 


^T" 


This  equation  can  be  combined  with  the  thermodynamic 
equation 


(13) 
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to  yield 

Taking  the  dot  product  between  fj V   and  equation  (2) 
leads  to 

JlV^  JTV-  VVa4  JLST  iSL*" + 

77  v  V<£  -  W-V7T  -O.  (15) 

Multiplying  (12)  byd>  and  utilizing  the  relationship 

rrv-v$>  =  y-(TTV4>)-  cJ)v-(ttv) 

results  in 

Substituting  the  identities 

d)^r  -  /i&si  -  T-^  Use 

and 

into  (15)  and  adding  (15)  to  (12)  yields  a  kinetic  energy 
equation  of  the  form 
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£Jx  v)t v- fe vav+  7n£ \y) +  J>_( it-  v3^  +  7rcj* 


dt 


c)^r        c)t  <>t  d^~ 


The  first  law  of  thermodynamics  applied  to  adiabatic 
conditions  is 

Cp  dr  -  d&£   -  o 
*t  at 

Substituting  p=  n  \J  +Pt  and  expanding  the  total  derivatives 
into  Eulerian  form  for  ^J~  coordinates  leads  to 

i__  (Cpf)  +  v- vCcpt)  +  ^-d^fiD  -  <x 

Multiplying  this  equation  by  yp  and  equation  (12)  by 
CpT,  then  adding  the  resulting  equations  gives  the 
potential  energy  equation 

1_(  FCpT)  -t-  v-  (rrCpT  vj+ AiirCpT  t)-  tt*Vwls) + 


^7T  7")    f- 


V- V(W) +sr<l£K2f]  •=.  o  . 


(17) 


•7- 


Combining  equations  (16)  and  (17)  gives 

dt  x  a 

^  a.  <^t  <)^r  (18) 

Thompson  ["4"]  outlines  the  following  three  main  aspects 
to  be  considered  in  solving  the  primitive  equations: 
formulation  of  boundary  conditions,  formulation  of  initial 
conditions  and  numerical  integration  of  finite  difference 
equations,, 

The  method  of  numerical  integration  utilized  in  this 
study  will  be  the  method  of  repeated  time  extrapolation 
outlined  by  Thompson  and  shown  to  be  computationally 
stable  for  a  simple  linear  model  with  a  grid  interval  the 
size  of  the  one  used  in  this  study  if  the  time  interval 
is  about  10  min0  or  less.   This  method  employs  a  centered 
difference  of  the  form 


<h  =  cb     +   (&&\    a  At. 


Initial  conditions  for  solution  of  the  primitive 
equations  present  a  real  problem.   For  meteorologically 
significant  large-scale  motion  there  is  a  near  balance 
between  the  coriolis  and  pressure  forces „   This  means  the 
resultant  acceleration  is  an  order  of  magnitude  smaller 
than  either  of  these  two  forces.   Small  errors  in  the 
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measurement  of  initial  wind  and  pressure  fields  may  lead 
to  intense  divergence  fields  whichj  in  turn,  give  large 
pressure  changes,  thus  creating  rapid  accelerations „ 
These  spurious  accelerations  create  fictitious  gravity- 
waves  o   To  exclude  these  waves  initially  in  this  study, 
a  form  of  non-divergent  initial  winds  is  used*   However, 
this  does  not  necessarily  exclude  gravity  waves  for  all 
time* 

The  forecast  region  used  in  this  study  covers  the 
major  portion  of  the  Northern  Hemisphere «,   Boundary 
conditions  for  a  limited  region  such  as  this  are  difficult 
to  define  since  the  proper  type  of  boundary  conditions 
depends  on  the  mathematical  character  of  the  primitive 
equations,  which  have  a  shifting  dual  nature,  sometimes 
elliptic  and  sometimes  hyperbolic .   In  this  paper  all 
parameters  are  assumed  constant  with  respect  to  time  on 
the  boundaries «   It  is  recognized  that  this  assumption 
may  lead  to  the  existence  of  spurious  development  near  the 
boundaries;  however,  it  is  felt  that  if  the  spatial  finite- 
difference  schemes  are  such  that  in  the  mean  potential  and 
kinetic  energies  are  conserved,  any  boundary  instability 
will  tend  to  be  averaged  out  over  the  forecast  region,,   This 
averaging  process  will  also  tend  to  suppress  spurious  gravity 
waves  resulting  from  truncation  and  round-off  errors . 

It  has  been  shown  by  application  (see,  for  example, 
Haltiner  [6])  of  Gauss'  equation  that  when  an  equation 


of  the  same  type  as  (1&)  is  integrated  over  a  closed 
volume  the  local  time  change  of  the  sum  of  kinetic  and 
potential  energy  vanishes.   It  follows  that  any  finite- 
difference  scheme  should  in  the  mean  conserve  the  sum  of 
the  potential  and  kinetic  energies.   Specifically,  this 
consistency  should  be  achieved  separately  with  respect  to 
both  vertical  and  horizontal  differencing. 

The  vertical  consistency  is  developed  as  follows. 
First  equation  (15)  is  applied  at  the  \T   level.   The 
vertical  differencing  is  then  performed,  leading  to 

dt-  *  AST     <*> 

ah-  j\-f-  '         Arr'Wt 


AT  dt 


AT' 


(19) 
AT 


(^,-4W=oJ 


since  ^  =0,  V^    =0  and  V2^V1V3 °      ^he  deriva^ive  -^ 
is  approximated  by  a  one-sided  difference.   A  similar 
application  of  equation  (17)  results  in 


^  +  (EjV,)j  =  o 


(20) 
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Adding  equations  (19)  and  (20)  and  comparing  the  result 
with  equation  (IS)   leads  to  the  conclusion  that  if  the 
consistency  discussed  in  the  previous  paragraph  is  to  be 
maintained,  certain  terms  of  equation  (19)  must  cancel 
corresponding  terms  of  equation  (20),   This  inference 
leads  to  the  following  equivalences: 


±  I*,  +  $f  *)]  . 


A  similar  procedure  at  the  ,\JT  level  utilizing  the 
identity 


v"(|£)  =    ^M  _d> 


leads   to 


tt  ^     ft  T3  =-  (f)a  +  d>5  -  acj)w 

and 

From  the  previous  equations  the  following  relationships  can 
be  obtained., 

T2=(i±fflk)  TX  +  lgi)kT  • 

tTh^F^-      2  1^3/    ■> 

Substituting  P  =200  mb,  p1=400  mb,  p  =£00  mb 
Xj  -*-  j 

k=0.2£6j/g°k,  T{  =h   and  C  =1.003 j/g°k  in  the  previous  system 
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of  equations  leads  to 

T2=  *643  TX  +  .410  T3  (21) 

<f>,  =  <£<+  +  .442  RTx  +  .530  RT3  (22) 

$2  =  -»05S  RTx  +  .530  RT3  (23) 

(p3=  4>h+   .053  RTX  +  .22  RT3  (24) 

If  these  relationships  between  geopotential  and  temperature 
are  maintained  then  the  vertical  differencing  scheme  will 
not  violate  the  energy  conservation  equations „ 

With  respect  to  simple  derivatives  of  the  form  ^jdl 
involved  in  (15),  the  centered  difference  approximation 

will  obviously  vanish  when  summed  over  interior  points  of 
the  grid  and  thus  would  in  the  mean  conserve  potential  and 
kinetic  energy.   However  this  does  not  take  care  of  the 
boundaries,  nor  terms  which  are  not  of  the  simple  form 
shown  above. 

Since  numerical  integration  is  to  be  performed  by 
repeated  time  extrapolation,  the  velocity  tendencies  (#)  - 
(11)  will  have  to  be  computed  from  q) ,  T  and  Va.  after  they 
have  been  extrapolated  forward  in  time.   This  will  require 
tendency  equations  for  7T~~and  T  since  (p  values  are  expressed 
in  terms  of  T  through  equations  (22)  -  (24) «. 
Integration  of  (12_)_  leads  to  the  approximate  form 

m  &  --Lv-Ljrfy-i-^O*  (25) 

at,  a 
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From  (12) 

substituting  for<^HT  from  (25)  gives 

^  i  ^[^(rrv^Trv,)!  • 


The  relationship 

utilized  in  equation  (17)  yields  for  the  ^  level 

ircp^r;  +  cPT;/Air^v-r/r^)/t-CpTT^-\7r  -t- 


(26) 


The  bracketed  terms  of  equation  (26)  may  be  replaced  by  the 
following  relationship 


Dividing  (26)  by  JT  Cp  gives 


(27) 


A  similar  approach  at  the  ^3  level  will  result  in  the 
following  tendency  equation 


£ll  +  V-  7T3  +  ^.(T68T3  -  T2 


03  3  T- 


Jt 


/r 


=0 


(23) 
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based  on  equations  (21)  -  (24)  with  <fu assumed  zero. 

These  tendency  equations  complete  the  system  required 
for  repeated  time  extrapolation*, 
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3.   Numerical  experiments  and  results „ 
1.   Energy  Considerations 

In  the  basic  model,  the  energy  balance  with  respect 
to  vertical  differencing  is  maintained  through  the 
relationship  between  temperature  and  geopotential  developed 
in  Section  2.   With  respect  to  horizontal  differencing, 
the  problem  may  be  considered  in  two  parts, 

A.  internal  energy  balance, 

B.  boundary  effects. 

For  the  former  problem,  partial  derivatives  of  a 
"closed"  form,  such  as^_£«  ,  can  be  maintained  in  balance 
through  use  of  a  simple  central  differencing  scheme.   As 
the  computation  proceeds  across  the  grid,  each  term  is 
added  and  subtracted  once  at  each  grid  point,  with  no 
change  in  total  energy.   However,  a  few  derivatives  with 
variable  coefficients  must  be  evaluated  in  the  horizontal 
acceleration  equations,  (3)  through  (11),  and  also  in  the 
temperature- change  and  stability-change  equations.   In 
these  terms,  of  the  form  v  .)£!,  the  coefficient  varies  as 
the  grid  is  traversed,  so  that  the  weighting  of  the  term 
at  a  specific  grid  point  changes  between  addition  and  sub- 
traction in  central  differencing.   This  can  lead  to  a  change 
in  total  energy  across  the  field,  with  resultant  build- 
up of  instability. 

The  following  approaches  were  tested  on  these  "open" 
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forms : 

a0   Terms  left  in  this  form* 

This  lead  to  a  comparatively  short  set  of  forms „   For 

example,  the  bracketed  portion  of  equation  (S)  was 

simplified  to  vl°ul  -  ul^  ul  -  v1f 0   However ,  this  form  did 

=}y     <>x 
not  fulfill  the  requirements  for  energy  balance ,  and  the 

model  became  unstable  in  a  very  short  time  with  overflow 

occuring  in  less  than  four  hours „ 

b»   Terms  placed  in  a  "closed"  form  as  a  sum  of  squares. 

The  bracketed  form  of  equation  (3)  now  become 

This  form  was  a  distinct  improvement  on  the  "open"  form; 
however,  variable  coefficients  involved  in"^,,  were  not 
affected .   In  any  case,  there  was  some  instability  and  this 
model  overflowed  in  about  six  hours „ 

Co   Open  form  with  multipliers  replaced  by  average 
values. 

In  this  model,  the  equations  of  the  first  approach 
were  used,  but  a  "central  average"  replaced  the  variable 
coefficient  in  each  case.   For  example ,  vg)u  became,  in 
finite  difference  notation: 

(v/.t-h)    +V,t:r-ft    J.  (u,(iTiQ   r.Utari)— J 
2  2^y 

This  model  gave  the  best  results  of  the  three,  and  was 

used  for  the  remaining  experiments „ 

The  problem  of  boundary  effects  proved  much  more 

difficult.   Due  to  computer  storage  limitations,  this  model 
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was  programmed  using  an  octagonal  grid  whose  southern 
boundaries  are  as  much  as  13°  from  the  equator „   Advection 
at  these  boundaries  leads  to  major  problems „ 

The  following  boundary  forms  were  tested^ 

a0   All  time  derivatives  set  equal  to  zero  along  the 
boundaries .   Since  the  boundary  values  were  held  constant , 
while  values  at  adjacent  points  were  allowed  to  vary, 
large  gradients  resulted,  giving  large-amplitude  gravity 
waves  which  moved  inward  from  all  sides  at  high  speed „ 
(See  Fig,  1) 

b„   Boundary  values  averaged  at  short  time  intervals. 

Several  smoothing  techniques  were  tried  which 

involved  only  the  outside  three  grid  points •   Since  a 

gravity  wave  cannot  cover  the  distance  between  two  grid 

points  in  less  than  one  hour  of  forecast  time,  this 

averaging  was  carried  out  once  each  hour0   The  most 

successful  technique  involved  placing  an  average  of  the 

three  outer  points  at  the  next-to-edge  point,  then  averaging 

this  value  with  that  at  the  outside  point  and  placing  the 

result  at  the  outside  point .   For  example; 

(Boundary)     A  B         C 

becomes       A+(A+B+C)/3         (A+B+C)/3      C 

2 

This  method  restricted  the  instability  due  to  boundary 

effects  to  the  outer  two  grid  points;  however,  in  spite  of 

the  heavy  averaging,  large  differences  developed  between 

the  outer  three  grid  points  after  about  five  hours  fore- 
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cast  time. 

c.   All  horizontal  velocities  set  equal  to  zero  along 
the  boundaries. 

Since  the  instability  developed  along  the  boundaries 
can  be  traced  to  energy  advected  into  or  out  of  the  forecast 
area,  zeroing  the  horizontal  winds  at  the  boundaries, 
thereby  zeroing  the  advection,  should  solve  the  problem . 
In  practice,  it  was  found  that  zeroing  the  horizontal  wind 
components  at  both  forecast  levels  at  the  outer  two  grid 
points  was  effective  in  controlling  boundary-produced  gravity 
waves.   Surprisingly,  this  did  not  lead  to  instability  in 
the  wind- component  fields  themselves  in  this  model .   This 
method  of  boundary  condition  control  was  used  throughout 
the  remaining  experiments „ 
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2.   Friction  Terras 

Two  friction  terms  were  developed  for  use  with  the 

horizontal  acceleration  equations  (£)  through  (11).   For 

horizontal  friction,  a  simple  Laplacian  multiplied  by  lCK£m/jec 

was  subtracted  at  each  time  step.   For  vertical  friction, 

the  terms  developed  by  Smagorinsky F  5j   were  used,  with  an 

average  skin  friction  coefficient  of  2.2  x  10~3  after 

Cressman  [  3~).   According  to  Smagorinsky,  the  vertical 

friction  at  any  level  may  be  computed  from 

vF  =  -ggir=  g_  (Tk+i  "X-i). 
K    dp       ap 

Assuming  7o=0, 

vF  =g&  ;    vF  =g__  (%-Tz) 

1  ^p  -?  AP 

Vertical  stress,  y   ,  is  taken  as  a  function  of  the  vertical 
wind  shear,  so  that 

or  since  K  is  proportional  to  the  magnitude  of  the  wind, 
VI  i  and  (•    is  a  constant  for  any  roughness  length 

i 


%  =(k  jv[v| 


^d  depends  on  roughness  length,  and  from  Cressman' s  paper 

2  -3 

an  average  value  of  2.2  x  10  J   was  used.   Level  4  was  taken 

to  be  the  top  of  the  friction  level  for  this  term,  with 

JV  ,  estimated  as  a  fraction  of  the  vectorially  extrapolated 

wind  by  the  equation 

Ivlf  =  0.36  (  (1.192u,-0.192u  )2+  (1.192v.-0.192v. )  )  . 
'  '  4  $  1  *  1 
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For\)^  the  cross  isobar  angle  at  level  4  was  taken  to  be 
a  function  of  latitude  alone  by  the  equation 

ctnj"=  l-10272f"  j    f  in  flec  • 
This  gives  a  cross-isobar  angle  ranging  from  near  zero 
at  the  pole  to  approximately  35°  at  the  base  of  the  grid 
(1$N),  which  is  in  close  agreement  with  average  values • 
The  magnitude  of  \/  was  obtained  from  the  geostrophic 
approximation  and  the  (pn    field,  which  in  turn  was  estimated 
by  a  simple  extrapolation  from  (h    and  $3  ,  namely 

$t   =1.192^-0.192^  , 
For  *Y ,    the  equation, 

Xm  -vk)2  gA  (V-v,) 

c  ^p    1  J 

can  be  readily  obtained  from  the  basic  T*  equation »   Here 
gf>  /flkP=l/az,  where  z=5.236  km,  the  standard  depth  of  the 
#00-  to  400-mb  layer.   For  (fK)2,  the  Rossby  and 
Montgomery  estimate  of  $00  gm/m  sec  was  used. 

This  series  of  equations  was  programmed  and  tested 
on  the  model.   However,  the  terms  involved  were  so  small 
compared  to  the  rest  of  the  model,  and  the  computations 
so  time-consuming  (approximately  50%  increase  in  computing 
time),  that  time  was  not  available  for  a  complete  evaluation. 
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3 .   Prognosis  of  Surface  Pressure 

The  prognosis  of  surface  pressure  was  one  of  the  most 
difficult  problems  encountered  because  the  surface  pressure 
change  of  the  basic  model  was  apparently  very  sensitive  to 
small  errors  in  velocity  aloft .   In  factj  the  surface 
pressure  field  served  as  a  diagnostic  chart ,  since  any 
instability  was  apparent  there  long  before  it  appeared 
aloft. 

Four  different  methods  were  tried  for  obtaining  the 
prognostic  surface  pressure  field: 

a.  Expansion  of  the  basic  equation  (25). 

The  resulting  pressure-change  field  was  extremely 
sensitive  to  changes  in  velocity  divergence  at  levels  1  and 
3,  and  large  pressure  changes  occurred  in  two  to  three  hours . 

b.  Averaging  around  the  grid  point c 

In  this  approach ,  the  products  of  the  pressure  value 
and  the  sum  of  the  wind  components  at  levels  one  and  three 
were  differenced  across  the  point  being  computed „   In 
finite  difference  form, 

TRjt.J  (V,G"i  ^(i+o)   -'"O'-O  faff-0  ^Q-ij  J  J  * 

This  model  was  not  quite  as  sensitive  to  divergence 
errors,  but  still  gave  unreasonable  surface  pressure  changes 
over  periods  as  short  as  6  hours . 

c.  Computation  of  surface  pressure  directly  from  the  cb 
fields . 
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In  this  case,  the  deviation  of  the  surface  pressure 
from  standard  was  estimated  at  each  time  step  from  the 
deviation  of  the  <p  fields  from  standard  at  300  and  400  mb 
and  the  hydrostatic  equation^   i.e.,  if  (p^PiP   is  defined  as 

K'-Qtmmo     '   etc'  then 

dp  2. 

and  if  we  let  th         -=.  rb         »  then  from  the  hydrostatic 
V\ooo        tsfc 

equation,  Api^p^^  and  p=pstci  +  Ap. 

This  model  was  programmed  and  computed.   However, 
since  the  (p    fields  are  based  on  temperature  and  the  basic 
inputs  to  this  model  were  the  temperature  patterns  at  400 
and  #00  mb,  the  resultant  surface  pressure  field  lacked 
some  of  the  usual  detail. 
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4  ,    Conclusions  and  suggestions  for  further  study- 
In  this  series  of  experiments ,  three  finite  difference 
forms,  three  boundary  forms,  and  three  surface  pressure 
forms  were  tested .   It  was  found  that  this  model  is 
extremely  sensitive  to  finite  difference  forms  used  and 
boundary  conditions  assumed  <,      However 5  when  a  new  form  gave 
a  closer  approximation  to  energy  balance ,  stability  was 
improved  in  every  case,  with  the  improvement  noted  being 
in  direct  proportion  to  the  approach  to  energy  balance,, 

It  is  the  opinion  of  the  writers  that  a  stable 
primitive-equation  model  of  this  type  can  be  achieved ,  and 
will  eventually  become  the  standard  form  of  numerical 
prognosis .   However,  faster  computers  with  much  larger 
storage  capacity  must  be  used  to  realize  this  goalo   At 
least  two  more  levels  of  temperature  input  are  needed  to 
delineate  the  geopotential  fields  clearly s  and  the  grid 
should  be  expanded  to  cover  at  least  a  whole  hemisphere 
to  reduce  boundary  problems 0   Even  with  high  speed  magnetic 
drum  storage ,  such  a  model  would  have  a  prohibitively  long 
running  time  on  a  computer  such  as  the  CDC  1604;  but  it 
might  be  completely  feasible  on  the  next  generation  of 
computers o 

For  the  present,  further  experiments  with  surface 
pressure  forms  are  suggested ,  possibly  using  some  simpli- 
fied version  of  the  equations  currently  under  development 
by  Arakawa  f 7J  .   Also,  the  friction  terms  could  be 
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streamlined  and  tested,  and  a  surface  terrain  term  included 
However,  any  major  revision  of  this  model  will  require  the 
use  of  a  larger,  faster  computer  than  any  currently 
available. 
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APPENDIX  I 

LIST  OF  FIGURES 
Figure 

1.  Section  of  6-hour  $00-mb  temperature  prognosis  with 
no  boundary  correction. 

2.  Section  of  12-hour  SOO-mb  temperature  prognosis, 
edge-averaged. 

3.  Section  of  l£-hour  SOO-mb  temperature  prognosis, 
boundary  advection  zeroed. 

4.  400-mb  Phi  analysis 

5.  400-mb  Phi  1$   hour  prognosis 

6.  400-mb  Temperature  analysis 

7.  400-mb  Temperature  IS  hour  prognosis 
3.   300-mb  Phi  analysis 

9.   SOO-mb  Phi  l£  hour  prognosis 

10.  800-mb  Temperature  analysis 

11.  £00-mb  Temperature  IS  hour  prognosis 
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0     00Z  20  JAN  94  PH  800 

PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60'  NORTH   LATITUOE 
SCALE    1:60,000.000  V. 


Chart  no.  6-B      1 


18HRPRQG  FROM, 
0    002  20  JfilN  6*i  Phi  800 

PROJECTION:  POUR  STEREOGRAPHIC-TRUE  AT  60"  NORTH  LATITUDE  . 
SCALE:   1:60,000.000  V- 


♦       18  HOUR 
880  T.PROdFRQM    ,  00Z'I120  JAN  64 

PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60"  NORTH  UTITUDE 
SCALE:   1:60.000.000 


0     00Z  20  JAM  64  PH  400 

PROJECTION     TOIAH  SIFREOGRAPHIC-TRUF  AT  60    NORTH   LATIllint 
<!CAI.F     1:60,000  000  >•: 


ritri    NUMERICAL  WEATHER   FACILITY 
MONTI  RFY    CAlirORNJA 


400       T       GNPL     .  00Z  20  JAN  64 

h       PROJECTION:  POI  AR  STEREOGRAPHIC-TRUE  AT  60    NORTH  LATITUDE 
9CALE    1:60.000,000  fcj*. 


;■:■♦■/.  18  HOUR 

400  T  PROG  FROM    .  00Z  -20  JAN  M 

PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60°  'NORTH  LATITUDE 
9CALE;    1:60.000,000  if. 
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